Magnetic Instabilities and Phase Diagram of the 
Double-Exchange Model in Infinite Dimensions 

R S Fishman*, F Popescu , G Alvarez*, J Moreno^, Th Maier*, and M Jarrell* 



* Materials Science and Technology Division, 
Oak Ridge National Laboratory, Oak Ridge, TN 37831-6032 



< 

\q . "Physics Department, Florida State University, Tallahassee, FL, 32306 



'Computer Science and Mathematics Division, 

Oak Ridge National Laboratory, Oak Ridge, TN 37831-6032 
up ! 

^Physics Department, University of North Dakota, Grand Forks, ND 58202-7129 and 

^Deparment of Physics, University of Cincinnati, Cincinnati, OH 45221 

s 

£ 

O 

> 

On 

m 
^f 

o 
\o 
o 

^-* 

B 

O 

o 
> 

X 



Abstract 

Dynamical mean-field theory is used to study the magnetic instabilities and phase diagram of 
the double-exchange model with Hund's coupling Jh > in infinite dimensions. In addition to 
ferromagnetic (FM) and antiferromagnetic (AF) phases, the DE model also supports a broad class 
of short-range ordered (SRO) states with extensive entropy and short-range magnetic order. For 
any site on the Bethe lattice, the correlation parameter q of a SRO state is given by the average 
q = (sin 2 (#j/2)), where 9i is the angle between any spin and its neighbors. Unlike the FM (q = 0) 
and AF (q = 1) transitions, the transition temperature of a SRO state with < q < 1 cannot be 
obtained from the magnetic susceptibility. But a solution of the coupled Green's functions in the 
weak-coupling limit indicates that a SRO state always has a higher transition temperature than 
the AF for all fillings p below 1 and even has a higher transition temperature than the FM for 
0.26 < p < 0.39. For 0.39 < p < 0.73, where both the FM and AF phases are unstable for small Jh, 
a SRO phase has a non-zero transition temperature except close to p = 0.5. As Jh increases, the 
SRO transition temperature eventually vanishes and the FM phase dominates the phase diagram. 
For small Jh, the T = phase diagram of the DE model is greatly simplified by the presence of 
the SRO phase. A SRO phase is found to have lower energy than either the FM or AF phases 
for 0.26 < p < 1. Phase separation disappears as Jh — ► but appears for any non-zero coupling. 
For fillings near p = 1, phase separation occurs between an AF with p = 1 and either a SRO or a 
FM phase. The stability of a SRO state at T = can be understood by examining the interacting 
density-of-states, which is gapped for any nonzero Jh in an AF but only when Jh exceeds a critical 
value in a SRO state. 

PACS numbers: 75.40. Cx, 75.47.Gk, 75.30.-m 



I. INTRODUCTION 

With the renewed interest in manganites [1[ and the revitalized study of dilute magnetic 
semiconductors 2], attention has once again focused on the double-exchange (DE) model 
and its variants. Applied to the manganites, the DE model is customarily formulated with 
one local moment per site and a Hund's coupling Jh taken to be much larger than the 
electron bandwidth W. In studies of dilute magnetic semiconductors, on the other hand, 
the local moments are sparse and the exchange coupling is comparable to W. While the full 
phase diagram of the DE model as a function of Jh and electron filling p has been studied by 
several authors J3, ]% |5|, |y] , the properties of the DE model for small coupling constant remain 
very much in doubt. In particular, there have been conflicting claims about the presence 
of phase separation (PS) p and incommensurate phases 3, |4j, |y] for small Jh, when the 
ordered FM and AFphases are magnetically frustrated by a RKKY-like interaction between 
the local moments |7|. 

This paper uses DMFT to evaluate the magnetic instabilities and T = phase diagram 
of the DE model. Developed in the late 1980's by Muller-Hartmann J8( and Metzner and 
Vollhardt [9|, DMFT exploits the fact that the self-energy becomes independent of momen- 
tum in infinite dimensions, where DMFT becomes formally exact. Even in three dimensions, 
DMFT is believed to capture the physics of correlated systems including the narrowing of 
electron bands and the Mott-Hubbard transition llCi . Although DMFT has been widely 
applied to the DE model |4laUll2lUWUUIlilldl, until now there has been no 
complete treatment of the phase instabilities and T = phase diagram of the DE model for 
arbitrary J H and p. 

We shall study a system with a bare semicircular density-of-states (DOS) given by 



N (u) = (8/nW 2 )Re J W 2 /4: — u 2 . In real space, this DOS belongs to an infinite- 

ri 

dimensional Bethe lattice or an infinite Cayley tree with no closed loops [19|. A finite- 
dimensional Bethe lattice with coordination number z c = 4 is sketched in Fig.l. Although 
the Bethe lattice lacks translational symmetry, it is quite convenient for calculations. As 
shown in Fig.l, the Bethe lattice can be partitioned into A and B sublattices so that both 
ferromagnetic (FM) and antiferromagnetic (AF) long-range order are possible. Due to the 
bounds ±W/2, the DOS of the Bethe lattice more closely resembles the DOS of two- and 
three-dimensional systems than does the unbounded DOS of the hypercubic lattice. Indeed, 




FIG. 1: A Bethe lattice with z c = 4 nearest neighbors. As shown, the Bethe lattice may be 
partitioned into A and B sublattices, denoted by the blue and red dots. 

pathological results have been obtained on a hypercubic lattice due to the asymptotic free- 
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As we shall see, the Bethe 



dom of the quasiparticles in the tails of the Gaussian DOS 

lattice also has the advantage that analytic results are possible in the limit of small Jh, 

precisely the regime where controversies persist. 

The high-temperature, non-magnetic (NM) phases of the Heisenberg and DE models have 
a correlation length £ that vanishes as z^— * oo. By contrast, the short-range ordered (SRO) 
states introduced in an earlier paper [20] possess some of the same characteristics as spin 
glasses: local magnetic order and exponentially decaying magnetic correlations J2l|. The 
SRO states are characterized by the correlation parameter q = (sin 2 (#j/2)), where 9{ are the 
angles between any spin and its neighbors. Overall, the neighboring spins describe a cone 
with angle 2 arcsin( v /g) around the central spin. This disordered state reduces to the FM and 
AF phases when q = and 1, respectively. Unlike the FM and AF phases, a SRO state with 
< q < 1 has extensive entropy and short-range but not long-range magnetic order. Since 
it is not possible to associate a wavevector with any SRO state, its transition temperatures 
TsRo(p,q) cannot be obtained from the magnetic susceptibility. Rather, T SRO (p, q) can be 
solved from coupled Green's function equations. Remarkably, a SRO state can be found 
with higher transition temperature and lower energy than the ordered FM and AF phases 
in a large region of phase space. 



The FM and AF transition temperatures may be obtained by solving either the coupled 
Green's function equations or the Bethe-Salpeter (BS) equations for the magnetic suscep- 
tibility The BS equations break up into separate relations for the charge and spin degrees 
of freedom. The vertex functions of the BS equations are consistent with the functional 



derivative of the self-energy with respect to the Green's function [17[. It follows that the 
DMFT treatment of the DE model is ^-derivable and hence thermodynamically consistent 
22j |. meaning that identical results are obtained from both ways of evaluating Tq and T$. 



This paper is divided into 7 main sections. In Section II, we formulate the BS equations 
for the susceptibility and derive the vertex functions. The thermodynamic consistency of 
the DMFT of the DE model is discussed in Section III. The BS equations are used in Section 
IV to derive implicit expressions for the FM and AF transition temperatures. Using coupled 
Green's functions, we rederive the transition temperatures of the ordered phases and evaluate 
the transition temperature of the SRO states. We also explicitly solve those relations in the 
weak-coupling limit of small Jh. A complete map of the phase instabilities of the DE model 
is provided in Section V. In Section VI, we present the T = phase diagram of the DE 
model both with and without the SRO states. Finally, Section VII contains a summary and 
conclusion. Appendix A provides results for the transition temperatures and T = energies 
of the SRO states in the weak-coupling limit. In Appendix B, we derive the self-consistency 
relation for the SRO states. Appendix C provides a convenient and general form for the 
kinetic energy of the SRO and ordered phases. 

II. BETHE-SALPETER EQUATIONS AND THE VERTEX FUNCTIONS 

The Hamiltonian of the DE model (also called the FM Kondo model when Jh < oo) is 
given by 

H = -tJ2 (c\ a c ja + c] a c ia ) - 2 J H 51 s, • Si (1) 

where c ia and Ci a are the creation and destruction operators for an electron with spin a at 
site i, Jh > is the Hund's coupling, Sj = cl a cr a pc i 0/2 is the electronic spin, and S; = Srrii is 
the spin of the local moment. Repeated spin indices are summed. The local moment will be 
treated classically, which is only a fair approximation for the magnanites with spin S = 3/2 
but a better approximation for dilute magnetic semiconductors with S = 5/2. For large z c , 



he hopping energy t scales as \j ^fz~ c and the local effective action on site may be written 

HI 

S cS (m) = -T^2c 0a (iv n ){G (iv n )~p + Ja af3 ■ m}c 0/ g(w n ), (2) 

n 

where J = JuS, v n = {2n + l)7rT, co a (iz/„) and co a (iv n ) are now anticommuting Grassman 
variables, m is the orientation of the local moment on site 0, and G_ (ii> n ) is the bare Green's 
function containing dynamical information about the hopping of electrons from other sites 
onto site and off again. 

Because S'eff(m) is quadratic in the Grassman variables, the full Green's function 
Q(m, iu n ) a p at site for a fixed m may be readily solved by integrating over the Grass- 
man variables, with the NM result 

g(m,iv n ) = \G Q [iv n ) l+Ja-m\ = — — — — , (3) 

LrQ\lV n ) — J 

where Go(w n )a/3 = 5apGo{iv n ) is diagonal above Tq. Averaging over the orientations of the 
local moment, (G(m,iv n ) a p) m = G{iv n ) a p = 8 a pG(iv n ) where 

G{t "' :) = G (»w)--7' - (4) 

If P(m) = Trf (exp(— S e s(ra)))/Z is the probability for the local moment to point in the 
m direction, then the average over m is given by (C(m)) m = j dQ m P(m)C(m). The 
effective partition function Z involves both a trace over the Fermion degrees of freedom and 
an average over m: Z = J dn m TiA exp(— S'eff(m)) J. Above Tq, -P(m) = l/An is constant. 
Consequently, the NM self-energy is Tj(iu n ) = G (iis n ) _1 — G(iz/ n ) _1 = J 2 Go(^n)- 

Likewise above Tq, the lattice Green's function Q(m, k, iv n ) a p has the expectation value 

(£(m, k, iv n ) a p) m = 5 a/3 (5) 

z n ~ e k - J 2 G {iv n ) 

wit h z„ = iv n + \i. Summing Eq.© over all k and equating the result to Eq.(JU), we obtain 



m 



m 



w 2 

Go{iv n ) = z n - —G(iu n ), (6) 



where W = ktyfz~ c is the full bandwidth of the semi-circular DOS 
Our goal is to evaluate the general susceptibility [15| 



E4£(^>n) = / dre^Y,^' {R ^ o) {( T rcl(T)c^T)clc os ) - (cU)^)}, 

l,n ° i L J 

(7) 
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FIG. 2: Diagrammatic representations of (a) the magnetic susceptibility XWkJ^ l {^ v li^ v n)i (b) the 
vertex function T^' K (ivi,iv n ), and (c) the Bethe-Salpeter equations. 



where u m = 2miiT and (C) involves both a Fermion trace and an average over all the 
rrij. Sketched in Figs. 1(a) and (b), the general susceptibility and the full vertex function 



-tT)fl,TU . 



iui,iu n ) are related by the lattice BS equation: 

/3»,<5k _ (0)fia,SK _|_ (OJ i/9a,^»7pw,ri/ vr 



^■C\,tUJ> 



±.<\,IU) 



■C\,%Un 



•,6k 



(8) 



Each term in this equation is a matrix in Matsubara space and matrix multiplication involves 
a summation over intermediate Matsubara frequencies, as shown in Fig. 2(c) [15|. The local 
susceptibility y^ a ^ K is obtained from the lattice susceptibility ^ a ' SK by summing over q and 
satisfies the local BS equation: 



/3a,5n _ (0)/3q,<5k , (0) fe^ftTi' yVT,5n 



lLO m C±.iu) 



(9) 



Since the vertex function is independent of momentum in infinite dimensions, it appears 
identically in both the lattice and local BS equations |2j,|2J|. 

For the bare susceptibilities, the electrons in each Green's function are uncorrelated, 
which means that each Green's function is separately averaged over m. So the bare, lattice 



susceptibility is given by 

T 

An. 
k 

The bare, local susceptibility is obtained by summing this expression over all q: 



xtl^ 5K ( iv u iv n ) = -jr^in Y.(G{m, k ' iu i)as)m(Q(rn, k + q, iv l+m ) K p) m - (10) 



XijJ? a ' SK (i u h iu n) = -T5i n G(iui) aS G(iiyi +m ) Kl3 , (11) 

where G{ivi) a p = 5 a pG{ivi) is the m-averaged, local Green's function given by Eq.(jlJ). It is 
also straightforward to evaluate the full, local susceptibility: 



Ba,6K / ■ ■ \ 

xLi \w h iv n ) 



T5 m ,o{(G(m,Wn)paG(m,ii>i)s K ) m - {G(m,iVn)p a }m(G(m,m)s K } m } 

-T8 ln (G(m, iui) aS G{m, m +m ) K p) m , (12) 



where Q{vn.,ivi) a p is the in-dependent, local Green's function given by Eq.(j3J). 

Due to the matrix structure of Eqs. (jlOHT2*|) . the charge and spin degrees of freedom 
decouple: a^x^- l3a ' 5K ^ a i K an d 5 al3 x^ /3a ' SK ^ 5 &K with analogous relations for the 

— q,2o; m - — q, 2aj m 

local susceptibilities. We define the charge and spin susceptibilities and vertex func- 
tions by X c ■ = 'Wx^'W^ X s ■ = < 3 X PaM <s/^ EL = S af3 T^6 KS , and 
EL m = a a8^-c^'ujl a K5- Similar definitions are employed for the local susceptibilities. Since 
Eqs. fllO|) and (|TT|) imply that the bare, charge and spin susceptibilities are identical, we shall 
dispense with the superscripts c and s on those quantities. 

After this decoupling, both the lattice and local BS equations separate into two sets of 
equations for the charge and spin degrees of freedom: 

X a ■ =X {0) +X (0) r?, X a ■ , (13) 



x a = y 



(0) +y (0) IX x a , (14) 



where a = c or s. The same charge and spin vertex functions enter the lattice and local 
equations. The spin BS equations with a = s were previously solved [15j in the strong- 
coupling limit J ^> W and J 3> T. 

To generally solve Eq.(JT3J) for the charge and spin lattice susceptibilities, we first obtain 
the charge and spin vertex functions from the local susceptibilities given by Eq. (jl4j) . The 
vertex function may be written as 

r £f K (^' ip n) = A ln af3 a ■ a 5 J mfi + B l m a 5a ■ <Tp K 8i n , (15) 



An = qn : ^ , , , (16) 

3(1 + C) kb n 

B i 13 J 2 Go(iui)Go(iui +m )aiai +m 

where a n = G Q (iv n )~ 2 — J 2 , b n = G (^ n )~ 2 — J 2 /3, and C = — (2J 2 /3) J2 n b^ 1 . In terms of 



A\ n and B m , the charge and spin vertex functions are then given by T^ m (iui,iiy n ) = 6B m S[ n 
and Y s iuJm {iv h iv n ) = AA ln 5 mfl - 2B l m S in . 

III. THERMODYNAMIC CONSISTENCY 

A thermodynamically consistent theory produces identical results from either the Green's 
function (which is on the single-particle level) or the partition function (which involves 
interactions on the two-particle level). Baym and Kadanoff .22] demonstrated that a 
sufficient condition for a theory to be thermodynamically consistent is that it be $- 
derivable, which means that a functional $({G(ir/ n )}), constructed from compact diagrams 
involving the full Green's functions and the bare vertex functions, satisfies the condition 
S(wn)a/3 = $Q/3G(iv n ) a p. The partition function of a ^-derivable theory satisfies the rela- 
tion 

- logZ = $ - £Tr f {£(^ n )G(^ n )} +Y,Trtlog{G(iu n )}, (18) 

n n 

which is stationary under variations of G_(iu n ). Whereas Baym and Kadanoff considered 
systems of interacting Fermions and Bosons, their ideas were later extended to systems of 
interacting electrons and spins [2a] and to disordered alloys [261 ] . The $-derivability of the 
DMFT of the DE model was first demonstrated in Ref . [l7j . 

The bare vertex function T\J^ a ' K (iui, iu n ) may be associated with the two-particle inter- 
action in the purely electronic effective action [27 1 

S' eS = -T^c 0a {iv n )G Q {iv n y l c 0a {iv n ) 

n 
T 3 

— j- Yl Safari + iuj m )c 0l 3(iVn)TfJf aM (iis h iv n )c$ K {ivi)c^(ivi + iu m ). (19) 

Ln,m 



Because the Fermion operators anticommute, the bare vertex function must satisfy the 
crossing symmetries r^_ ( {ivi,ii'i +m ) = r>J ( _ n {iu n+m ,iv n ) = -T)J m [%v u %v n ). There 
are two ways to evaluate the bare vertex function. First, we can take the J — ► limit of 



the full vertex function rf^ re given by Eq. (jl5|) . Alternatively, we can associate the lowest- 
order, J 2 contribution to the partition function Z = J <ifi m Tr f (exp(— S^m)) j with the J 2 
contribution to the partition function Z' = Tr f (exp(— S' cS ) J. Both methods yield the same 
result: 

Liu^'^fah Wn) = o^^ 2 {°"/3a ' a sJm,0 ~ O'Sa ' (TpJln], (20) 

which does indeed satisfy the required crossing symmetries. 

But S' eS produces an inequivalent infinite-dimensional theory than produced by 5 , e ff(m), 
as seen in the different expansions of Z and Z' to order J 4 : 



Z = Z \l - J^Go^n) 2 + ^J 4 E^o(^) 2 Go(^n) 2 + 0(J 



(21) 



Z' = Z (l-J 2 EGoK) 2 + ^ 4 EGoN 2 GoK) 2 + 0(J 6 )). (22) 

Therefore, the DE action cannot be replaced by a purely electronic action containing only 
two-body interactions and it is not possible to construct a diagrammatic expansion in powers 

ofr<°>. 

Nevertheless, we can construct a functional <&({G(iv n )}) satisfying the condition 
S(ii/ n ) oa = 5$/5G(iv n ) aa . Starting from Eq.(j2J and Dyson's equation for the self-energy, 
we obtain 5Y,(ii>i) aa /5G(ii> n )pp = (K' 1 )^ + 5i n 5 a pG(iv n )~ 2 , where 

Kg = J^P^U = -sM 2 -(- + KSJ + -^(28 a , - 1). (23) 

Inverting this Jacobian yields the general result 

gN» x j2a l\ 2 na„\A j2 aian fox A 

6G{iv n )pp 3b n [2a n -36 n J 3(1 + C) kb n \ > 

= -TTtZ%(m^n), (24) 

where the right-hand side is the full vertex function of Eq. (jl5j) . The functional $ must 
exist because the vertex function is well-defined and the curl of the self-energy vanishes: 
5Y,{iv n ) aa / 5G{ivi) pp — 5T,(ii/i) pp / 'SG(iu n ) aa = 0. So despite the absence of a diagrammatic 
perturbation theory, the DMFT of the DE model is ^-derivable and hence, thermodynami- 
cally consistent. 



10 



IV. TRANSITION TEMPERATURES 

Since the electronic spin-susceptibility is given by x ol (q) = Yli,nXq,iu m =o^' Ul ^ Un )' the 
condition that x^Cq) diverges is equivalent to the condition that Det(x s . _ n ) _1 = 0. To 
obtain the FM or AF transition temperatures, we require the q = or Q susceptibilities with 
uj m = 0. The AF wavevector Q is defined so that the Fermi surface is nested at half filling 
or €k+Q = — £k- Although the Bethe lattice sketched in Fig.l lacks translational symmetry 
and nonzero wavevectors q are ill-defined, the AF wavevector Q may be associated with the 
bipartite nature of the Bethe lattice. To be safe, the Neel temperature will be evaluated in 
two different ways, only one of which relies on the existence of an AF wavevector. 

For either FM or AF ordering, the electron filling p can be related to the chemical potential 

/i through the relation 

32 

p = 1 + T£ReG(^ n ) = 1 - — T^Rei? n , (25) 

n V" n 

where Goiii^n) 1 = z n + R n and the complex function R n satisfies the implicit relation 

W 2 
a n R n + —(z n + R n )=0. (26) 

lo 

Since a n = (z n + R n ) 2 — J 2 , R n falls off like — W 2 /(16z n ) for large n. Notice that p — 1 when 
\x = 0. Since all of our results are particle- hole symmetric about p = 1, we shall restrict 
consideration to the range < p < 1. 

As discussed in Section II, the BS equations decouple into charge and spin relations. 
Since the charge vertex function T^ (ivi, iv n ) = 6B l m 5[ n is diagonal in Matsubara space, the 
condition Det(x c ■ _ n ) _1 = has no solution for q = or Q. In other words, there is no 
charge- density wave instability within the DE model. 



A. Ferromagnetic ordering 

When q = 0, the bare, zero-frequency lattice susceptibility is given by 
(o) ,. . , 4T( zi - J 2 G Q {iv{) \ 

w V \](zi - J 2 G {tvi)) 2 - W 74/ 
With the help of Eq. (j26|) . we find that the full spin-susceptibility at q = satisfies 

- -^T^W,^ + *)'<* + ^ ~ J ^> + «'/ 3 )} - w^c]w,' (28) 

11 



where bi = (zi + Ri) 2 — J 2 /3. Consequently, the condition that Det(x s _ _ ) x = can 
be written as an implicit condition for the Curie temperature: 

3 ^(z n + R n y(z n + 2R n )-J2(z n + 4R n /3) ' l ; 



which must be supplemented by Eq. (|25|) for the electron filling. 

If J/W > 1/4, the interacting DOS splits into two bands [lq. In the lower band centered 
at —J with < p < 1, all electrons have their spins parallel to the local moments; in the 
upper band centered at J with 1 < p < 2, all holes have their spins parallel to the local 
moments. As J/W increases, both upper and lower bands are narrowed by correlations and 
their widths approach W = W/y2. To recover the large J limit, the chemical potential 
must be written as p, = Jsgn(p — 1) + 5p where —W/2 < S/i < W/2. Eq. (f29| then reduces 

to y,iij 

R 2 

? rJ - 31V732 = lj (30) 



where R n is now given explicitly by R n = (Jz 2 — W 2 /8 — z n )/2. So in this limit, T c is 
proportional to W and is independent of J. Due to the strong Hund's coupling, T c vanishes 
at half filling (p = 1) or when the lower band is completely full because electrons are unable 
to hop to neighboring sites. 

In the weak-coupling regime of T ^C J ^C W, the summation in Eq. (|29jl may be replaced 
by an integral over frequency u, as shown in Appendix A. Carefully treating the branch cut 
along the real axis between ±W/2, we find that 

T c = H^l - (2^/W) 2 {(2^/W) 2 - 1/4}, (31) 

which is proportional to J 2 /W as expected in the RKKY limit. For small filling (// close to 
— W/2) and a fixed J, Tq initially increases with p and then begins to decrease, vanishing at 
\i = W/A corresponding to a filling of p ~ 0.39. Because Tc(p, J) is particle-hole symmetric, 
it is the same for p electrons (2 — p holes) and 2 — p electrons (p holes) per site. These 



results are in qualitative agreement with Nolting et al. 28, 29], who evaluated the Curie 



iemperature of the s — f model with quantum local spins using both equations-of-motion 



I and modified RKKY J29J techniques. Nolting et al. obtained FM solutions for small 
J in a concentration range significantly narrower than found here, probably because their 
approach retains the non-local fluctuations produced by the RKKY interaction. 



12 



For small fillings and e = 2fi/W + 1 < 1, T C W/J 2 « 8^2£ 1 / 2 /(37r). Since p « 
8A/2e 3 ' /2 /(37r), T c grows like p 1 / 3 for small p. 

B. Antiferromagnetic ordering 

This section presents two calculations for the AF Neel temperature: one based on the 
magnetic susceptibility and the other on coupled Green's functions. Their equivalence was 
established in Section III, which demonstrated that the DMFT of the DE model is thermo- 
dynamically consistent. While the first technique explicitly assumes the existence of an AF 
wavevector Q, the second does not. 

1. From the susceptibility 

Evaluating the Neel temperature from the susceptibility requires the bare, zero-frequency 
lattice susceptibility for q = Q: 

x^kuo =-w*[ 1 *, - j*gom — r n - (32) 

After some algebra, we find that the full spin-susceptibility at q = Q satisfies 

= -2(35 ln - %—{ Zl (z l + Ri) 2 - J 2 {z t + 2^/3)j - ^ °'°" . (33) 

Consequently, T N is found from the condition that Det(x* . _ ) _1 = or 

-—Y ^ = 1 (34) 

3 ^z n (z n + R n ) 2 -J 2 (z n + 2R n /3) ' l ; 

which bares some similarity to Eq. (j29j) for Tc- 
In the strong-coupling limit, Eg. (J34J) reduces to 

R? 

? Rl + 3lV 2 /32 = *' (35) 

with a plus sign multiplying 3W 2 /32 compared to the minus sign in Eq. (j3()j) . It is straightfor- 
ward to show that Tn < for all p in this limit. Hence, the AF phase instabilities disappear 

as J/W — > oo. 

13 



In the weak-coupling regime, we may once again replace the Matsubara sum by an inte- 
gral, with the result 



8J 2 { /l + \A-(2/i/W0 2 \ , 4 . 3/2 ^ 

T » = MH 2|H/V ) - \^^W- 3(1 - (2,/Wr) }. (36) 

For p < 0.73, Tn is negative so the system becomes FM for < p < 0.39. 

Close to half filling, T N given by Eq.()3b|) becomes positive and the DE model supports 
an AF phase. We find that Tn > for 0.73 < p < 1, in which range Tq given by Eq.(|31|) is 
negative. The divergence of Tn in Eq. (J36|) as p — > 1 or \x — > signals a breakdown in the 
weak-coupling expansion and the appearance of a gap in the interacting DOS, as discussed 
later. 

2. From coupled Green's functions 

The previous sub-section used the nesting condition for the AF wavevector Q whereas, 
strictly speaking, nonzero wavevectors are not defined on the Bethe lattice. Hence, it is 
worthwhile to re-evaluate Tn using a Green's function approach that skirts the question of 
whether an AF wavevector exists. This Green's function approach will also come in handy 
when considering the SRO state later in Section IV. 

Dividing the Bethe lattice into A and B sub-lattices, we obtain coupled relations between 
the bare and full local Green's functions: 

fiftuO -1 = z n I_ - ^£M(i„ n ), (37) 

lb 

where 77 =| or { on the A or B sites and 77 is opposite to r\. If the bare inverse Green's 
function is parameterized as 

d$\iv n Y x = (z n + R n )L + Q^gi z , (38) 

then Eq. (|37|) implies that R n and Q^ are solved from the condition 

RnL+Q^a, = -^{{(z n + R n )l+a- (jm + Q^z)} _1 )_. (39) 

To linear order in the sublattice magnetization, R n and Q^ satisfy the implicit relations 

W 2 
a n R n = — — (z n + Rn), (40) 
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^=^{^( i+ £) +jj ^}- (41) 

where (m)^ = M^z = ±Mz, the sign depending on whether r\ =\ (A site) or 77 =| (B site). 
The linear relations for Qffl may be readily solved, with the result 

Vra (z n + i?„) 2 z„-J 2 (^ n + 2i? n /3)- l J 

After integrating exp(— S^m)) over the Grassman variables, we find that the probability 
for the local moment on sublattice rj to point along m is 



P(m)W oc expj^logf 1 - 2J Q* m A I K exp(pj eff M M m z ). 



(43) 



The last relation defines the effective interaction: 



J ff = - 2TJ V ^- = -2TJ 2 V — ('441 

which is the same on each sublattice. Finally, Tn is solved from the condition M^ = 
J e ffM^f3/3, which reproduces Eq.(|34|). So the earlier results for T N do not depend on the 
existence of an AF wavevector! 



C. Other magnetic phases 



Earlier studies of the Heisenberg model on the Bethe lattice [30|, |3l| proposed a spiral 
magnetic state with a wavevector that characterizes the twist of the spin between nearest 
neighbors. But it is quite easy to show that a Bethe lattice can not be partitioned into more 
than two sublattices. For example, if we try to partition a Bethe lattice with z c = 4 into 
three sublattices, then sites C and A (red and purple dots in Fig. 3 (a)) would be next-nearest 
neighbors but the two sites C (purple dots) (presumably on the same magnetic sublattice) 
are actually next-nearest neighbors of each other. So unless the central site is singled out, a 
Bethe lattice does not support a magnetic state with more than two sublattices. 

More recently, Gruber et al. [32j argued that this restriction may be removed as z c —> oo 
because the paths with the largest weight always move outward from the central site. This 
analysis uses the fact that the surface of the Bethe lattice contains a non- vanishing fraction 
of its sites [33]. Since the period-three state constructed in Ref. 32j can only be formed 
around a central site, it is not isotropic. Such states will not be considered further in this 
paper. 
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q = 3/4 



FIG. 3: (a) A possible partitioning of the Bethe lattice into three sublattices with blue, red, and 
purple sites, (b) A colinear SRO state on a z c = 4 Bethe lattice with q = 3/4, meaning that 1/4 
of the neighbors of any site have the same spin while 3/4 have the opposite spin. The two spins 
are denoted by blue and red dots. 



The broad class of SRO states first introduced in an earlier paper 20] is characterized 
by the correlation parameter q = (sin 2 (#j/2)), where the average is taken over all neighbors 
of any site on the Bethe lattice. Since the azimuthal angles fa of the neighboring spins are 
assumed to be random with (exp(ifa)) = 0, the neighboring spins will on average sweep out 
a cone with angle 2arcsin( Ai /g) about the central spin. Due to the topology of the Bethe 
lattice, a SRO phase cannot be ordered (except in the FM or AF limits when q = or 1). 
For the special case of a colinear SRO state, 0, = or i\ on every site and q is the fraction 
of neighbors of any site that have the opposite spin. An example of a colinear SRO state 
with z c = 4 and q = 3/4 is sketched in Fig. 3(b). 

Quite generally, the SRO phase has short-range but not long-range magnetic order. If M n 
is the average magnetization of the lattice sites at a distance na (lattice constant a) from a 
central site, then it is easy to show that M n = (1 — 2g)M n _i. So the magnetization decays 
exponentially like \M n \ = \M \ exp(— na/!;) with the correlation length £ = —a/log \2q— 1| > 
0. The correlation length diverges in the FM and AF limits and vanishes in the NM state 
with q = 1/2, when the spins around any site are completely disordered. A SRO phase can 
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have either FM (0 < q < 1/2) or AF (1/2 < q < 1) correlations, as discussed further below. 

Also working on the Bethe lattice, Chattopadhyay et al. [5] assumed that all of the 
angles 9i were identical and characterized the resulting state by its "incommensurate spin 
correlations." We avoid this description because nonzero wavevectors different are ill-defined 
on the Bethe lattice and the X(q) formalism of Muller-Hartmann J8[ cannot be adapted 
from the hypercubic to the Bethe lattice. In our view, the SRO phase is more accurately 
characterized by the magnetic correlations around any site. 

In the absence of the SRO states, a large portion of the magnetic phase diagram is NM. 
This poses a problem because Nernst's theorem requires that the entropy vanishes at T = 0. 
Generally, the entropy of a SRO state is difficult to evaluate. But in the z c — > oo limit, the 
entropy of a colinear SRO phase is 



•J col 



-q\ogq-(l-q)log(l-q), (45) 



which vanishes only in the FM or AF phases. Of course, the entropy of a non-colinear SRO 
phase must be greater than S co \. Therefore, the formal dilemma posed by Nernst's theorem 
is not resolved by the presence of SRO states. As discussed in the conclusion, entropic 
ground states frequently arise in infinite dimensions. In finite dimensions, the SRO state 
may evolve into a state with incommensurate SRO. 

Since a SRO phase with correlation parameter q between and 1 has no long-range 
order, the transition temperature Tsro(p, q) cannot be obtained from the condition that the 
magnetic susceptibility diverges. But we can easily solve for T^ro {Pi q) by generalizing the 
Green's function technique used in sub-section IV. B. 2. With the SRO state defined above, 
the self-consistency relation between the bare and full local Green's functions is derived in 
Appendix B: 

G^(^)- 1 = z n I- ^{qQW(iv n ) + (1 - q)GW(w n )}, (46) 

where G} v ' and Gy' are spin-reversed Green's functions. The quantization axis t] is defined 
by the spin on the central site 0. The above relation was previously obtained in Ref. [5J] 
under the more restrictive condition that d\ is constant. Once again parameterizing the bare 
inverse Green's function by Eq. ljHHj) . we find that R n is independent of the magnetization 
M = (m z ) on site (assumed spin up) and is given by Eq- fjlUj) while Qffl = —Qx' has the 
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FIG. 4: The correlation parameter q (solid) and the associated maximum SRO transition tem- 
perature Tq^q {p) (dashed) versus filling p for the SRO states in the weak-coupling limit. Also 
plotted in the thin short-dashed line are the weak-coupling limits of Tq for q = and Tn for q = 1. 



linearized solution 



(2q - l)MJa n R % 



(47) 



(z n + R n f(z n + 2(1 - q)R n ) - P(z n + 2(2 - q)R n /3) ' 
which reduces to Eq . (|42j) when q = 1. 

The effective interaction J e yj is again given by the first equality in Eq. (j4"4"|) so that 
Tsro{p, q) is solved from the implicit relation 



2J 2 



(2g-i)E 



Rn 



(48) 



3 " ' '' ^ (z n + 2(1 - q)R n )(z n + R n f-J 2 (z n + 2(2- q)R n /?>) 
which reduces to both the FM and AF results, Eq. (f2H]) and (pMJ) . when q = and q = 1. It 
also follows that TsroG 5 ? 9) = in the NM state with g = 1/2. 

In the strong-coupling limit W <C J and T <t^ J, Tsro(p, q) is solved from the condition 



E 



^n 



(49) 



i22+3W*7(32(2<z-l)) 

It is straightforward to show that T SRO (p, q) is negative for 1/2 < q < 1 and is positive but 
always smaller than Tc(p) for < q < 1/2. So as expected, the SRO states are unstable in 
the strong-coupling limit. 
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Results for Tsro(p, q) in the weak-coupling limit are given in Appendix A. Maximizing 
Tsro(p, q) with respect to q, we obtain Tg^o (p) plotted in Fig.4. Remarkably, a SRO phase 
with AF correlations can always be found with a higher transition temperature than the AF 
state. As q — *■ 1, the AF and SRO transition temperatures meet at p = 1 or half filling. Even 
in the range of fillings between about 0.26 and 0.39, where Tq is nonzero, a SRO phase with 
FM correlations or < q < 1/2 has the higher critical temperature! Also notice that the 
SRO transition temperature is nonzero in the range of fillings between 0.39 and 0.73, where 
neither the FM nor AF states was stable for small J/W. The point at which q = 1/2 and 
Tsro (p) = li es slightly below p = 0.5. Recall that q changes discontinuously at T^f (p) 
from 1/2 in the NM phase above to a value greater than or smaller than 1/2 in the SRO 
phase below. 

V. MAGNETIC INSTABILITIES OF THE DOUBLE-EXCHANGE MODEL 

We shall first present results for the magnetic transition temperatures in the absence of 
the SRO states. These results are then revised to include the SRO states. Bare in mind 
that just because a FM, AF, or SRO state has the highest transition temperature for a given 
J/W and p does not imply that it remains stable down to zero temperature or even that it 
is reached at all: a pure state may be bypassed by the formation of a PS mixture. The issue 
of PS can only be resolved by calculating the free energy as a function of filling. We shall 
perform such a study for T = in Section VI. 

In Fig. 5, the Curie and Neel temperatures are plotted versus filling for a variety of cou- 
plings. When p > 0.39, a critical value of J is required for Tq to be nonzero. So the 
weak-coupling regime is restricted to electron concentrations smaller than 0.39. As J/W in- 
creases, the AF instabilities are restricted to a narrower range of concentrations near p — 1. 
For large J/W at half filling, the nearest-neighbor hopping can be treated perturbatively 
and superexchange leads to a Neel temperature proportional to W 2 / J. So the sliver of AF 
phase stabilized around half filling rapidly narrows with increasing J/W. By contrast the 
FM phase dominates as J/W increases. In the limit J/W — > oo, the Neel temperature 
vanishes even at p = 1. 

For a fixed filling, Tq and T-^ are plotted versus J/W in Fig.6. As expected, the weak- 
coupling curves track the numeric results for small J/W but deviate for larger J/W. Whereas 
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FIG. 5: The Curie and Neel temperatures versus filling p for a variety of couplings J/W. 



0.015 
0.010 
0.005 



■p=0.05 

• p=0.05 (wc) 

p=0.3 

p=0.3 (wc) 



(a) 



0.1 



0.03 


- p=0.85 


p=0.95 ( b ) 
p=0.95 (wc) 






0.02 


,-''' """' 


~"~ 


0.01 


.'/' .<^""'" 




n 


i.^-r" i . i . i . i . 



0.2 0.3 

J/W 



0.4 



0.5 



FIG. 6: The (a) Curie and (b) Neel temperatures versus J/W for several fillings. Also shown are 
the weak-coupling transition temperatures obtained using perturbation theory. 

the Curie temperature reaches a plateau at a value of order W, the Neel temperature for 
any p < 1 vanishes at a finite value of J/W. Again, the filling range of AF instabilities 
rapidly narrows as J/W increases. 

The phase instabilities of the DE model in the absence of SRO states are plotted in 
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FIG. 7: The phase instabilities of the DE model (a) assuming only FM, AF, and NM phases and 
(b) also including the SRO phase. 

Fig. 7(a). The NM region shrinks as J/W increases and is absent for J/W > 0.5. Since the 
interacting DOS splits for J/W > 1/4 .16], the NM region survives even after the formation 
of sub-bands. Whereas the AF region persists at half filling for all J/W, the sliver of AF 
phase rapidly narrows with increasing coupling. The FM/NM phase boundary in Fig. 7(a) 



li 



agrees quite well with that found by Auslender et al. 

Now including the SRO states, we plot the correlation parameter q along with T^q (p) 
in Fig.8. For J/W > 0, a transition into the AF state (q = 1) occurs in a narrow range 
of fillings close to p — 1. This range vanishes both as J/W — > and as J/W — > oo. It 
has a maximal extent (from p = 0.78 to 1) when J/W ~ 0.33. For J/W < 0.33, Tgao (p) 
vanishes at a single filling that corresponds to q — 0.5, meaning that there is no SRO. For 
0.5 > J/W > 0.33, TgR£> (p) i s nonzero for all fillings. As can be seen in Fig. 8(a) for 
J/W = 0.375, the correlation parameter for 0.5 > J/W > 0.33 jumps from q = in the FM 
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FIG. 8: (a) The correlation parameter q, (b) the maximum transition temperature T<^q(p)/W, 
and (c) the normalized transition temperature Tq r q ,(p)W/J 2 for various values of J/W versus 
filling p. 

phase to a nonzero value 0.5 < q < 1 in the SRO phase. For J/W > 0.5, one of the ordered 
phases always has a higher transition temperatures than the SRO states and q jumps from 
to 1 at the FM/AF boundary. The deviation between the numeric and weak-coupling results 
in Fig. 8(c) as p — > 1 is quite expected: in the limit p — > 1, the weak-coupling result for the 
SRO transition temperature diverges, thereby invalidating the condition T <C J under which 
it was derived. 

The magnetic instabilities of the DE model including SRO states are mapped in Fig. 7(b). 
Notice that the region formerly labeled NM in Fig. 7(a) is now occupied by SRO states. The 
SRO phase also occupies the formerly AF region at very small J/W. As anticipated above, 
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the NM curve (where T^q(p) = and q = 1/2) intersects the FM phase boundary at 
J/W ~ 0.33. This curve separates a SRO region with AF correlations (1 > q > 1/2) on the 
left from a SRO region with FM correlations (1/2 > q > 0) on the right. 

VI. MAGNETIC PHASE DIAGRAM OF THE DOUBLE-EXCHANGE MODEL 

The presence of a continuous range of SRO states with correlation parameter q obviously 
complicates the magnetic phase diagram of the DE model and the issue of PS. Because the 



aflflQ, 



weak-coupling behavior of the DE model has been the subject of some debate 
we shall treat this limit analytically. We then provide numerical results for the interacting 
DOS that explain the stability of the SRO phase. Finally, we present the T = phase 
diagram of the DE model, both with and without the SRO states. 

A. Weak-coupling limit 

In their numerical work at T = 0, Chattopadhyay et al. |5( obtained a very complex 
phase diagram, with a range of SRO states above J/W ~ 1/16 and PS between AF and FM 
phases for 0.1 < p < 1 below this value. For small J/W, a pure FM phase was found to be 
stable for < p < 0.1. We now show that PS diappears in the exactly solvable limit of zero 
temperature and weak coupling. 

At zero temperature, the FM, AF, and SRO energies are evaluated after once again 

parameterizing Gq (z = iv + /i)" 1 by Eq. (|58j) . Then, R(z) is obtained from the quartic 

equation 

R{z + R)Uz + 2(1 - q)R) 2 - J 2 } + ^-(z + 2(1 - q)R) 2 = 0, (50) 

1 > lo 

which differs from Eq. (j26j) (except when q = 1/2) because we are now working in the broken- 
symmetry phase. Similarly, Q^\z) = —Q^'(z) is given by 

Q mU) = (2a-l) MJR{Z) (51) 

Q [z) [Zq ^ z + 2{l _ q)R{z y W 

which is likewise valid for any J. 

The potential energy can generally be written as (V)/N = —2JM(sq z ), where the expec- 
tation value of the electronic spin on site for any temperature is 

(s 0z ) = ± J dvteG«HzU*U = 7rI y 2(2 8 g _l) /^ g(T)( " ) - (52) 
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For any J, the potential energy may be written as 

1 ,„/ ^ !6^ 2 f , R{z) 

Of course, the NM potential energy is obtained by setting q = 1/2. 

Second-order results for (V(p,q))/N are provided in Appendix A. For small J/W, (s 0z ) 
is of order J/W so the electronic spins are not frozen at T = 0. Bearing in mind that 
Tsko(p, Q = 0) = Tc(p) and Tsro(p, q = 1) = T^(p), the difference between the potential 
energies of a FM, AF, or SRO phase and a NM is given by A(V(p,q))/N = —3Tsro(p,<i), 
where both sides are evaluated to order J 2 /W. So the potential energy is lowered when the 
transition temperature is positive. 

A general expression for the kinetic energy of the SRO phase is proven in Appendix C. 
At T = 0, Eq.dUl) becomes 

Jf( K ) = -^Y,J d "G M (zU{(l - q)G^(z) 0a + qG^(z) 0a }. (54) 

xi contrast to alternative expressions for the kinetic energy of the FM and AF phases 

lOJ . Eq.([54|) does not assume any particular wavevector for the ground state. Independent 

of q, the zeroth-order kinetic energy is given by E (p)/N = — 2W/(3tt)(1 — <5 2 ) 3//2 , where 

5 = 2fi /W and /i is the chemical potential evaluated to zeroth-order in J/W. For any J, 

the SRO kinetic energy is given by 

±-(K(p,q)) = -^- jdvR{z) 2 \l + (1 ~ 2g)j2 A (55) 

and the NM kinetic is recovered by setting q = 1/2. It is important to realize that the 
second-order kinetic energy contains two contributions: one from the dependence of the 
integrand on J/W for a fixed \x\ the other from the dependence of the chemical potential 
n(p) on J/W for a fixed p. 

Using the condition that p — 1 — 32/ (nW 2 ) f dvIleR(z) depends only on the filling and 
not on J, we may expand \i in powers of J/W as \i = fi + (J/W) 2 fi2 + ■ ■ ■■ Then fj,o(p) is 
solved from the condition 

p = l + - J6Vl - S 2 + sin -1 s\ (56) 

and H2{p) is given in terms of //o(p) by \ii = W5/{(1 — q) 2 + 5 2 {2q — 1)}. The divergence of 
//2 as q — > 1 and p — > 1 (or S — > 0) signals the opening of a gap in the interacting, AF DOS 
as discussed further below. 
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Incorporating the change in chemical potential, the result for the second-order, in- 
dependent correction to (K(p,q))/N is given in Appendix A. To second order in J/W, 
the difference between the kinetic energies of a FM, AF, or SRO phase and a NM is 
A(K(p,q))/N = 3Tsro(p> <z)/2- Hence, kinetic energy is lost (the electrons become more 
confined) when the transition temperature is positive. 

Consequently, the net change in energy E(p, q) = (K(p, q) + V(p, q)) in the weak-coupling 
limit is given by AE(p,q)/N = — 3Tsro(p, <?)/2, which includes the FM (q = 0) and AF 
(q = 1) phases as special cases. The same relation would be obtained for a Heisenberg 
model H = —(1/2) J2ij JijS% • Sj with classical spins within Weiss mean- field theory (which 
becomes exact in infinite dimensions), so that AE/N = — J(k = 0)S 2 /2 = — 3Tc/2. For 
large J/W, however, the DE value of AE(p, q = 0)/N in the FM phase is proportional to J 
whereas Tq is proportional to W. So the proportionality between AE(p, q)/N and Tsro{p, q) 
(and the analogy with the Heisenberg model) only holds in the weak-coupling limit. 

In the limit of small J/W, the ground state energy is minimized by the same correlation 
parameter q that maximizes the transition temperature! This result is not surprising: for 
small J/W, the transition temperature is also small so the correlation parameter that ap- 
pears at TsroG ? q) a l so minimizes the energy at T = 0. As a consistency check, the results in 
Appendix A may be used to show that the derivative of the energy is related to the chemical 
potential by dE(p, q)/dp = Np,, where both sides are expanded to second order in J/W. 

B. Interacting Density-of-States 

To better interpret the weak-coupling results, we construct the interacting DOS per spin 

"<<".«> = — I>c*"(* - » + *>» = Ji Im {«}„„«■ < 57 > 

As J/W — > 0, N(u, q) reduces to the bare DOS N (to) for any q. The AF DOS can be solved 
analytically for all J: 



^ 8U„ W 2 /4 + J 2 -lo 2 



which vanishes for \u\ < J and \u\ > JW 2 /A + J 2 . Hence, N{u,q = 1) contains a square- 
root singularity on either side of a gap with magnitude 2 J. As J/W — > oo, the width of each 
side-band narrows like W 2 /8 J. It is also possible to explicitly evaluate the DOS at u = for 
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FIG. 9: The interacting DOS N(ui,q) (normalized by 1/W) versus u/W for J/W = 0.2 in the 
FM or q = (solid, blue), SRO phase with q = 1/4 (long dash, light blue), NM or q = 1/2 (dash, 
green), SRO phase with q = 3/4 (small dash, red) and AF or q = 1 (dot dash, purple) phases. 



anyg: N(u = 0, q) = (8/irW 2 (l -q)) Re yV 2 - J 2 , which vanishes for J > J e = W(l-q)/2. 
'.■Pot a NM, the critical value of J required to split the DOS is J c = W/4, as found earlier 
lfij . For an AF at T = 0, J c = since a gap forms for any nonzero coupling constant. For 
a SRO phase, J c > 0. 

Taking J/W = 0.2, we numerically solve for R(z) and plot the interacting DOS versus 
u for q = 0, 1/4, 1/2,3/4, and 1 in Fig.9. The FM DOS has kinks at u = ±(W/2 - J). 
For energies between the kinks, both up- and down-spin states appear; on either side of 
the kinks, the bands are fully spin-polarized. These kinks disappear in a SRO state with 
q > due to the absence of long-range magnetic order. The width of the DOS shrinks as 
q increases from to 1. These results agree with the second-order change to the chemical 
potential given above. In fact, the divergence of the AF result for /i 2 as p — > 1 is caused by 
the appearance of a gap in the AF DOS for any nonzero J. 

For any temperature, coupling constant, and correlation parameter, the total energy of 
the DE model may be written as an integral over the interacting DOS: 



^E(p,q) = 2jdujLuf(u)N(u,q), 



(59) 



2(3 



where f(u>) = l/(exp(/3(uj — ji)) + 1) is the Fermi function and /i = nip). At T = 0, this 
relation may be expressed in the compact form 

><^4M 1+ ™)' (60) 

which agrees with Eqs.(|o^|) and (fo3j) . 

From these results, it is not hard to understand why a SRO phase can have a lower energy 
than an AF. Due to the narrowing of the AF DOS, the AF energy may be higher than that 
of a SRO state with J < J c (q) and no energy gap. If the chemical potential of the AF 
phase lies sufficiently far from the energy gap, then the broader DOS of a SRO phase with 
J < J c will favor that state over an AF. As p — > 1, this condition is impossible to satisfy for 
arbitrarily small J and the AF becomes stable. 

For a given p < 1, the expansions of the kinetic and potential energies in powers of ( J/W) 2 
are valid when J < J c (q) so that the DOS of a SRO state with correlation parameter q is 
not gapped. Conversely, for a fixed J, the expansions fail sufficiently close to p = 1 and 
q = 1 that a gap develops in N(u,q). 

C. Ground-state phase diagram and phase separation 

When only FM, AF, and NM states are considered, PS occurs due to the discontinu- 
ity in the chemical potential at the FM/NM and NM/AF phase boundaries. Since these 
discontinuities become weaker as J/W decreases, PS disappears as J/W — > 0. 

Once the SRO states are introduced, the physical energy E^ mm \p) is obtained by mini- 
mizing E(p, q) with respect to q for a fixed p. Although the chemical potential is a continuous 
function of filling, PS still occurs for non-zero J due to the formation of an energy gap in 
the AF state, so that (1 / N)dE^ mm ' / dp\ p= i- = —J. A necessary condition for PS between 
fillings p\ < 1 and P2 = 1 is that the second derivative of the energy .E( mm )(p) must change 
sign for some filling p* between p\ and p2- But if E^ mm \p) is expanded in powers of (J/W) 2 , 
then the condition that d 2 E^ 11 ^ (p) / dp 2 \ p * = may be written 

7T J 2 



(61) 



sVT^ W 2 5 2 

where 5 = 2hq/W and no is the chemical potential evaluated to zeroth-order in J/W. Since 
the value q* of q at p* is less than 1, the expansion in powers of (J/W) 2 is valid for J < J c (q*) , 
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FIG. 10: The T = energy E^ mm ' (p) and the value of q that minimizes the energy versus filling for 
J/W = 0.5. A Maxwell construction (given by the dots and dashed line) between fillings p\ = 0.647 
and p2 = 1 demonstrates the existence of PS. 

as discussed previously. When J/W — > 0, the right-hand side of Eq. lfrTTj) vanishes so the 
condition for PS cannot be met in the limit of vanishing J/W except as 5 — > or p* — ► 1. 
Near p — 1, 5 ~ (1 — p)ir/4 so Eq. (j61j) also implies that the difference p 2 — Pi grows linearly 
with J/W. These results contradict Chattopadhyay et al. (5|, who seem to use a Maxwell 
construction involving the difference E^ mm \p) — E (p) rather than the total minimum energy 

To show this more explicitly, we have used £ ,( - mm ) (p) to find the PS region through Maxwell 
constructions between fillings p\ < 1 and p 2 — 1. One such Maxwell construction is pictured 
in Fig. 10, where a phase mixture along the dashed line provides a lower energy than the 
pure phase on the solid curve. For J/W = 0.5, the SRO states are bypassed and PS occurs 
between a FM with p\ = 0.647 and an AF with p 2 = 1. 

The T = phase diagram of Fig. 11 (a) assumes that only the NM, FM, and AF phases 
are allowed. Without the SRO state, the ground-state phase diagram is fairly complicated. 
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FIG. 11: The T = magnetic phase diagram of the DE model (a) assuming only FM, AF, and 
NM phases and (b) also including the SRO phase. 

As implied by the argument above, PS disappears in the limit of vanishing J/W but the 
PS regions grow with increasing J/W. Two PS regions appear around the fillings of the 
AF/NM and NM/FM phase boundaries for J/W < 1. When < J/W < 0.16, the lower 
PS region in Fig. 11 (a) mixes an AF with p < 1 and a NM. For 0.16 < J/W < 0.28, an AF 
with p < 1 is unstable and the lower PS region mixes an AF with p = 1 and a NM. For 
J/W > 0.28, PS occurs between a FM and an AF with p — 1. 

Surprisingly, the T = phase diagram of Fig. 11(b) becomes much simpler once the 
SRO states are included. For J/W < 0.33, a SRO state is stable over a range of fillings 
that diminishes with increasing J. Near half-filling, PS is found between a SRO state with 
< q < 1 and an AF with p = 1. In agreement with the above discussion, the PS region 
grows linearly with J/W. For J/W > 0.33, the SRO phase is skipped and PS occurs between 
a FM with q = and an AF with p = 1. The special value J/W ~ 0.33 agrees with the 
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coupling where the NM line intersects the FM boundary in the instability map of Fig. 7(b). It 
also corresponds to the coupling where the AF instability covers the widest range of fillings. 

The PS region in Figs. 11 (a) and (b) diminishes with increasing J/W. In the strong- 
coupling regime, our results are in qualitative agreement with other authors J3j,[rj|. But for 
J/W = 2, we obtain PS over a somewhat larger range of fillings than obtained in Ref . jl3| . 
Aside from the behavior at very small J/W, the phase diagram of Fig. 11(b) agrees rather 
well with that presented by Chattopadhyay et al. p. In particular, those authors also found 
that the SRO phase disappears when J/W exceeds about 0.3. 

Together, Figs. 7 and 11 provide a bird's eye view of the magnetic instabilities and phase 
evolution in the DE model. In Figs. 7(a) and 11(a), a triangle of phase space remains NM 
for all temperatures down to T = 0; in Figs. 7(b) and 11(b), a somewhat larger triangle of 
phase space remains in the SRO phase down to T = 0. Keep in mind, however, that the 
ordered phases of Fi g. 7 may be bypassed when {J/W,p} lies in the PS region of Fig. 11. For 
example, Furukawa []j| found that for J/W = 2, an AF with p < 1 was excluded by PS 
between an AF with p — 1 and a NM. 

VII. CONCLUSIONS 

We have used DMFT to evaluate the magnetic phase instabilities and ground-state phase 
diagram of the DE model. Surprisingly, a SRO phase may have a higher transition temper- 
ature than the long-range ordered AF and FM phases. The ground-state phase diagram of 
the DE model actually simplifies in the presence of the SRO states. 

At first sight, the stability of disordered NM and SRO ground states may be hard to 
understand. After all, the presence of a disordered state at T = contradicts Nernst's 
theorem, which requires that the entropy of the local moments be quenched at zero tem- 
perature. However, entropic ground states appear quite frequently in infinite-dimensional 
calculations. Within the DMFT of the Hubbard model, AF ordering on a hypercubic lattice 
at half filling can be frustrated by hopping between sites on the same sub-lattice |l0j]. This 
additional hopping produces a NM Mott insulator at T = for large U/W, again in viola- 



tion of Nernst's theorem. Also on the hypercubic lattice, the DE mode 
state extending all the way to p = when J/W is sufficiently small 4j. 



has a NM ground 



Although they 



did not study such states in detail, Nagai et al. |4J detected the signature of incommen- 
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surate states for weak coupling. Therefore, many-body theories in infinite dimensions are 
frequently unable to lower the symmetries of disordered ground states. Of course, DMFT 
cannot access phases described by a non-local order parameter. It should also be noted that 
Nernst's theorem only applies rigorously to quantum systems and not to a DE model with 
classical local spins. 

For small J/W, the energy of the entropic SRO state is lowered below the energies of 
the FM and AF phases by the RKKY interaction |7|] , which generates competing FM and 
AF interactions [lj|. As the dimension is lowered, the SRO state may evolve into the state 
with incommensurate correlations (IC) obtained in one- and two-dimensional Monte-Carlo 
simulations by Yunoki et al. J3j. Replacing "SRO" by "IC" the phase diagram of Fig. 11(b) 
bears a striking resemblance to the phase diagrams of Ref.|3j. In agreement with our results, 
Yunoki et al. found that the short-range correlations in the IC phase change from AF to 
FM as p decreases. The phase diagram of Fig. 11(b) also looks quite similar to the three- 
dimensional phase diagram obtained by Yin p , who compared the energies of the AF and 
FM phases with the energy of a two-dimensional spiral phase carrying momenta along the 
(111) direction. Hence, the infinite-dimensional phase diagram of Fig. 11(b) is closely related 
to the results of finite-dimensional calculations. 



On a translationally-invariant hypercubic lattice [10| in infinite dimensions, the SRO 
phase may map onto an incommensurate spin-density wave (SDW) with a definite wavevec- 
tor. In that case, the transition temperature of the SDW could be obtained from either 
the coupled Green's functions or the magnetic susceptibility. However, there are no indica- 
tions for SDW ordering in the limited studies of the DE model on an infinite-dimensional 
nypercnbic lattice ,4, and we consider S nch a po S8 i bill t y to be remote. Mnch mote likely, 
the SRO phase would map onto a disordered state with short-range magnetic order and 
broad maxima at wavevectors q associated with the nesting of the Fermi surface, as found 
in Monte-Carlo simulations of one- and two-dimensional hypercubic lattices |3j. 

Bear in mind that DMFT studies of the DE model on a hypercubic lattice sacrifice a 
key advantage of the Bethe lattice: the possibility of obtaining analytic results due to the 
relationship, Eq.([S|). between the bare and full Green's function. For example, both sides of 
the expression AE(p, q)/N = — 3Tg RO (p, <?)/2 were evaluated independently to order J 2 /W. 
This gives us great confidence in our results for the energy and critical temperature in the 
weak-coupling limit. Such analytic results are not possible on a hypercubic lattice. Indeed, 
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pathological results have been found |16j] on a hypercubic lattice due to the asymptotic 
freedom of the quasiparticles in the tails of the Gaussian DOS. In most instances, more 
physical results are obtained on a Bethe lattice, where the bounded DOS more closely 
resembles the DOS in two and three dimensions. 

Working on a Bethe lattice has allowed us to clarify the behavior of the DE model at 
weak coupling. We have demonstrated both analytically and numerically that PS disappears 
as J/W — > 0. Whereas the FM phase is stable for p below about 0.26, a SRO state is 
stable above. As p — > 1, the SRO phase evolves into the long-range ordered AF phase. 
The transition temperature of the SRO phase exceeds that of the FM and AF phases for 
0.26 <p< 1. 

Whether the SRO phase is a new kind of spin glass or a spin liquid can only be resolved 
by future studies of the static and dynamic susceptibilities. An intriguing possibility is 
that the local moments are frozen within a spin-glass state but that the electronic spins 
are subject to quantum fluctuations within a spin-liquid state. Two of the most important 
unresolved questions about spin glasses are whether there exists a true thermodynamic spin- 

Hss transition 24| an d whether a model without quenched disorder can support a spin glass 
. We have answered those questions for the DE model on an infinite-dimensional Bethe 
lattice. Above Tsro(p, <?), the NM state is characterized by the absence of SRO around 
any site. Below Tsro(p, q), short-range magnetic order develops with a correlation length 
£(g) = — a/log|2g — 1| and the entropy, while still extensive, is reduced. Although not 
marked by a divergent susceptibility, the transition at Tsno(p,q) is characterized by the 
development of short-range magnetic order and a reduction in the entropy. Future work 
may establish connections between the SRO phase and the spin-glass solution to the random 



Ising model on a Bethe lattice 36j], where the Edwards- Anderson order parameter can be 
explicitly constructed. 

The strength of DMFT applied to DE-like models with a Kondo coupling 2Jrrij • Sj 
between the local-moment spin S* = 5nij and the electron spin Sj is that the Fermion 
degrees of freedom in the local action may be integrated exactly, leaving a function only of 
the local moment direction rrij. Consequently, formally exact results are possible for a range 
of models where the local moments are treated classically. This stands in stark contrast to 
the case for the Hubbard model [lOj], where the action cannot be analytically integrated and 
quantum Monte Carlo must be used to solve the local problem. We should also mention 
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that the classical treatment of the local spins is more questionable in the manganites, where 
S = 3/2, than in Mn-doped GaAs, where S = 5/2. 

The discovery of a SRO solution to the DE model in infinite dimensions has wide rami- 
fications. Our results open the door for future studies of the disordered ground-state of the 
DE model unencumbered by the numerical complexities of Monte-Carlo simulations [3J] in 
finite dimensions. 
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APPENDIX A: WEAK-COUPLING RESULTS 

This appendix presents the weak-coupling results for the SRO transition temperature 
and energy. For low temperatures, Y^ndi^n) — > (l/2vrT) / du g(y) so that 

J 2 (2g-1) fj R^v) 

v + R{v)f{v + 2(1 - q)R(v)) -J 2 (v + 2(2 - q)R(v)/3) ' 

(Al) 



Tsko{p, q) = — J dv ■ 



where v — iu + jjl. In the weak-coupling limit T <^i J <^i W, this becomes 

m , , J 2 (2g-1) r , R(z) /A , 

T SK o(p, q) = ^ / du f i m ^7 on ^TTV A2 

3vr J (z + R{z)) (z + 2(1- q)R(z)) 

where z = iv + jjlq and R(z) falls off like W 2 /z for large \z\. The zeroth-order chemical 
potential /io is defined in terms of the filling p by Eq. (|56|) . Using complex analysis and 
carefully treating the branch cut between ±W/2 along the real axis, we find that 

T SRo{M )W_ 32 l l dyy IT^y 2 -^-^l\ (A3) 



J 2 37r J 8 V • y 2 + b 

with 5 = 2/io/IU and b = (1 — q) 2 /{2q — 1). Clearly, Tsro(p, q) vanishes as q — > 1/2 or 

b — > oo. 
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In the regime < q < 1/2, the integral for Tsro{p, g) with b < gives 
T S Ro(p,q)W 8 



3g 



2^/l^Tq 



J 2 
Tan" 1 



97r(l-2g) 



1 3VT^^ - 4(1 - 5 2 ) 3/2 (l - 2q) 



[l-q)\6\-y/T=?q 



Tan 



-l 



l-g^l + v/T^V 



V^T^¥q J V VT^q 

This result can be extended to the regime 1/2 < q < 1 by using the identity 

J(l- q )\8\-^T=Tq' 



(A4) 



1 



Tan 



VT^Pq 

9 

log 



Tan 



-l 



'l-q)\6\+^T=2q-' 



\/T^2g 



VW^ 



q + VT^^/W^T 



5 2 {2q-l) 



qf 



(A5) 
(A6) 



These expressions yield Eq. (|?TT]) for Tq and Eq. (f3T)|) for Tn as q — > and g — > 1, respectively. 
Using Eg. (|54)l . it is straightforward to obtain the kinetic energy in the weak-coupling 
limit: 

2J 2 



^ K ^ = h Eo(p) + .W(l-2q) 



2VT^8 2 + 



Q 

jl^Tq 



Tan 



_ x {l-q)\8\-^T=Tq 



Tan 



-i 



;i- g )H + v /r^2f 



(A7) 



Of course, the NM kinetic energy is obtained by taking the limit q — » 1/2, with the re- 
sult (K(p,q = 1/2)) /N = E (p)/N + (1QJ 2 /3ttW)(1 - 5 2 f> 2 . The difference between 
the g-dependent and NM kinetic energies in the weak-coupling limit is then given by 
A(K(p,q))/N = 3T SRO (p,q)/2. 

Likewise, it is straightforward to evaluate the potential energy to order ( J/W) 2 , with the 
result: 

AJ 2 



1 
N 



(V(p,q)) 



*W{l-2q)\ V jT=2q- 



Tan 



J{l- q )\8\-y/T=?q 



Tan -1 



q)\5\ + ^T=Tq 



(A8) 



VT^Pq J V VT^q 

In the limit q^ 1/2, the NM result is (V(p, q = 1/2)) /N = -(32 J 2 /3nW)(l-5 2 ) 3 / 2 . Hence, 
the difference between the g-dependent and NM potential energies in the weak-coupling limit 
is given by A(V(p, q))/N = — 3Tsro(p, g)- It follows that the change in total energy is given 
by AE(p, q)/N = — 3Tsro{p, qO/2, which includes the FM (g = 0) and AF (g = 1) phases as 
special cases. 
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APPENDIX B: SELF-CONSISTENCY RELATION FOR THE SRO STATES 

Following the cavity method of Ref . [lOj , the self-consistency relation on the Bethe lattice 
between the bare and full Green's functions may be written 



/ 



G^lVn)" 1 = Z n l-t 2 Y.G{lV n )u, (Bl) 

i 

where t 2 = W 2 /16z c and G_{iv n )u is the local Green's function on site i. The sum on i is 
restricted to nearest neighbors of site 0. If the spin on site is ordered along the z axis, 
then we can write the local Green's function on site as 

G{iv n ) m =\ 9ln ° ] . (B2) 

V g in ) 

The creation and destruction operators on site i are taken to refer to a quantization axis 
that is rotated by angle 9{ (rotation axis in the xy plane at an angle <pi from the x axis) 
with respect to z. Then G(iu n )u = R(9i, <pi) G_{iVn)ooE.{9ii (piY is given by 

„,. , ( cos 2 (9 l /2)g 1n + sm 2 (9 l /2)g in isin0 ie -^(^ n -^ n )/2 \ 
G_\lVn)ii = ■ (B3) 

^ -isin^e^(^ Tn - g in )/2 cos 2 (6 i /2)g ln + sin 2 (^/2)# Tn J 

Summing over all z c neighbors of site and assuming that <pi is random so that 
(1/^c) T,' i exp(±i(p i ) = 0, we find 

i£Q-,,)„ = (l-4 9T " M^f 91 " M. (B4) 

z ' i \ gm) \ g,„) 

where q = (l/z c )Y^'iSin 2 (6i/2). Upon noting that the two Green's function matrices above 
are spin-reversed, we obtain Eq. 



APPENDIX C: KINETIC ENERGY 

Previous forms for the kinetic energy utilized the wavevector q = or Q of the FM or 
AF phase [l0|. In this appendix, we develop an expression for the kinetic energy that also 
applies to the SRO states, where the wavevector is not defined. Begin with the relation 

1 {K) = -^E/ n ^^Tr f ((cU a + cL Coa )e- 5 ), (CI) 



iV ZZjf 
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where Z$ = J IL; <if2 m .Trf (exp(— S)j and the sum over i is restricted to nearest neighbors of 
site 0. The action S = AS + So + S^ is separated into terms that involve site and the 

with site removed. The kinetic part of the action containing 
site is 

AS = -t J* dr J2 (ct a (r)c ja (r) + c] a {r)c Qa {r)) . (C2) 

Performing perturbation theory in the term involving the (Oi) bond and using t = WjkJ~z c , 
we find 

^<*> = ^E//^( c L(0)c ja (r))(cL(r)c 0a (0)) 

c i 

= ^ T E G ' ( ' 7) (^)^{(1 - q)G^\iPi)p a + qGMiivi)*,}, (C3) 

l,T) 

where the local averages at sites and i are performed using S e g of Eq.(j2J). Removing the 
single (Oi) bond from AS does not affect those local averages in the z c — > cxd limit. 
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